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The knowledge of the redshift of multiple images in cluster-lenses allows to determine precisely 
the total projected mass within the Einstein radius. The observation of various multiple im- 
ages in a same cluster is opening new possibilities to constrain the curvature of the universe. 
Indeed, although the influence of f2 m and £l\ on the images formation is of the second order, 
observations of many multiple images at different redshifts formed by a regular cluster-lens 
should allow to constrain very accurately the mass distribution of the cluster and to start to be 
sensitive to the cosmological parameters entering the diameter angular distances. We present, 
analytical expressions and numerical simulations that allow us to compute the expected error 
bars on the cosmological parameters provided an HST/WFPC2 resolution image and spectro- 
scopic redshifts for the multiple images. Numerical tests on simulated data confirm the rather 
small uncertainties we could obtain this way for the two popular cosmological world models: 
tt m = 0.3±0.24, n x = 0.7±0.5 or fi m = l.±0.33, Q x = 0.±1.2. Our method can be applied 
from now on, on real clusters. 



1 Introduction 

Recent works on constraining the cosmological parameters using the CMB and the high redshift 
supernovae seem to converge to a new "standard cosmological model" favouring a flat universe 
with Q m ~ 0.3 and £l\ ~ 0.7: White 1 and references therein. However these results are still 
uncertain and depend on some physical assumptions, so the flat £l m = 1 model is still possible 
(Le Dour et a/.H). It is therefore important to explore other independent techniques to constrain 
these cosmological parameters. 

In cluster gravitational lensing, the existence of multiple images - with known redshifts - 
given by the same source allows to calibrate in an absolute way the total cluster mass deduced 
from the lens model. The great improvement in the mass modeling of cluster-lenses that includes 
the cluster galaxies halos (Kneib et alM, Natarajan & KneibH) leads to the hope that clusters 



can also be used to constrain the geometry of the Universe, through the ratio of angular size 
distances, which only depends on the redshifts of the lens and the sources, and on the cosmolog- 
ical parameters. The observations of cluster-lenses containing large number of multiple images 
lead Link & Piercei (hereafter LP98) to investigate this expectation. They considered a simple 
cluster potential and on-axis sources, so that images appear as Einstein rings. The ratio of such 
rings is then independent of the cluster potential and depends only on Q m and assuming 
known redshifts for the sources. According to them, this would allow marginal discrimination 
between extreme cosmological cases. But real gravitational lens systems are more complex con- 
cerning not only the potential but also off-axis positions of sources. They conclude that this 
method is ill-suited for application to real systems. 

We have re-analyzed this problem building up on the modeling technique developed by us. 
As demonstrated below, we reach a rather different conclusion showing that it is possible to 
constrain £l m and £l\ using the positions of multiple images at different redshifts and some 
physically motivated lens models. 

Troughout this paper we have assumed Hq = 65 km s _1 Mpc -1 , however the proposed 
method is independant of the value of Hq. 

2 Influence of Q m and il\ on the images formation 



2.1 Angular size distances ratio term 
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In the lens equation: 9s 
contained in the term F 
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, the dependence on O m and Q\ is solely 
For a given lens plane, F(z s ) increases rapidly 
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up to a certain redshift and then stalls, with significant differences for various values of the 
cosmological parameters (see Fig. |I]). Thus in order to constrain the actual shape of F(z s ) 
several families of multiple images are needed, ideally with their redshifts regularly distributed 
in F(z s ) to maximize the range in the F variation. 
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Figure 1: Left. F(z s ) for zi 



0.3 and various cosmological models. Right. F(z S 2)/F(z a i) 
and for zi = 0.3, z 3 i = 0.7 and z 3 2 = 2. 



as a function of f2„ 



If we consider fixed redshifts for both the lens and the sources, at least 2 multiple images are 
needed to derive cosmological constraints. In that case F has only an influence on the modulus 
of 9\ — 9s- So taking the ratio of two different F terms provides the intrinsic dependence 
on cosmological scenarios, independently of Hq. A typical configuration leads to the Fig. |l| 
plot. The discrepancy between the different cosmological parameters is not very large, less than 
3% between an EdS model and a flat low matter density one. The figure also illustrates the 



expected degeneracy of the method, also confirmed by weak lensing analyzes, with a continuous 
distribution of background sources (e.g. Lombardi & Bertini ). 



2.2 Relative influence of the different parameters 

We now look at the relative influence of the different parameters, including the lens parameters, 

to derive expected error bars on Q m and Q\. To model the potential we choose the mass density 

distribution proposed by Hjorth & Kneiblil, characterized by a core radius, a, and a cut-off radius 

s S> a. We can then get the expression of the deviation angle modulus Dq 2 =|| 9\ — 9s \\- 

For 2 families of multiple images, the relevant quantity becomes the ratio of 2 deviation 

angles for 2 images 0n and 9i2 belonging to 2 different families at redshifts z s \ and z S 2- Let's 
D e 

define R$ n e I2 = J1 ■ With several families, the problem is highly constrained because a single 
potential must reproduce the whole set of images. In practice we calculate — ° n ' g/2 versus 
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the different parameters it depends on. We chose a typical configuration to get a numerical 
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evaluation of the errors on the cosmological parameters: zi = 0.3, z s \ = 0.7, z S 2 = 2, - — = 2, 

0/1 
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= 10 (9 a = a/DoL,9 s = s/Dol) and we assume O m = 0.3 and Q\ = 0.7. We then obtain 
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the following orders of magnitudes for the different contributions : 
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As expected, even with 2 families of multiple images the influence of the cosmological param- 
eters is of the second order. The precise value of the redshifts is quite fundamental, therefore 
a spectroscopic determination (dz = 0.001) is essential. The position of the (flux-weighted) 
centers of the images are also important. With HST observations we assume d9j = 0.1". 

So even if the problem is less dependent on the core and cut-off radii (in other word the 
mass profile), they will represent the main sources of error. Taking d9 a /9 a = d9 s /9 s = 20 %, we 
then derive the errors dVt m and d$l\ from the above relation in the flat low matter density we 
chose. We did this computation for different sets of cosmological models. Indeed the errors we 
will obtain with this method change significantly with respect to VL m and Q\ . All other things 
being equal apart from the cosmological parameters, we plot dQ m and dQ\ for a continuous set 
of universe models (Fig. ||). For instance in the 2 popular cosmological scenarios, we have : 

n m = 0.3 ±0.24 S! A = 0.7±0.5 or ft m = 1 ± 0.33 O A = 0±1.2 

As this can be easily understood from the Fig. [I] degeneracy plot, the method is in general 
far more sensitive to the matter density than to the cosmological constant, for which the error 
bars are larger. 

However the results we could obtain this way are as precise as the ones given by other 
constraints. But these errors are just typical; provided spectroscopic and HST observations, 
they depend mostly on the particular cluster and the potential model chosen to describe it. 
They could be quite tightened with a precise model, and by increasing the number of clusters 
with multiple images. 

3 Constraint on (O^, Oa) from strong lensing 

3.1 Method and algorithm for numerical simulations 

We consider basically the potential introduced in section 2.2. After considering the lens equation, 
fixing arbitrary values (f2^,f2°) and a cluster lens redshift zi, our code can determine the images 




of a source galaxy at a redshift z s . Then taking as single observables these sets of images as 
well as the different redshifts, we can recover some parameters (the more important ones being 
Co, 9 a or 9 S ) of the potential we left free for each point of a grid (£l m ,£l\). The likelihood of 
the result is obtained via a ^-minimization, where the \ 2 is computed in the source plane as 
follows : 

i=l j=l a Si 

The subscript i refers to the families and the subscript j to the images of a family. There is 
a total of Ya=i n l = N images. 9 l S j is the source found for the image 9\^ in the inversion. 9 l SG is 
the barycenter of the 9 l S j (belonging to a same family). Finally if is the error on the position 
of the center of 9\- and Aj the amplification for this image, then a^ Si = cP u j \J~^-. 

3.2 Numerical simulations in a typical configuration 

To recover the parameters of the potential (i.e. do, 9 a , 9 S and adjusted lens parameters), we 
generated 3 families of images with regularly distributed source redshifts. 

For starting values {Sl° m , tl° x ) = (0.3, 0.7) we obtained confidence levels shown in Fig. ||. The 
method puts forward a good constraint, better on Q, m than on Q\, and the degeneracy is the 
expected one (see Fig. Q). Concerning the free parameters, we also recovered in a rather good 
way the potential, the variations being Ado ~ 150 km/s, A9 a ~ 3" and A9 S ~ 20". 

This is an "ideal" case, of course, because we tried to recover the same type of potential we 
used to generate the images, the morphology of the cluster being quite regular and the redshift 
range of the sources being wide enough to check each part of the F curve. Such a simple approach 
can be applied to regular clusters like MS2137-23, which shows at least 3 families of multiple 
images including a radial one, see Fig. But the spectroscopic redshifts are still currently 
missing for this cluster. 

4 Conclusions & prospects 

Following the work of LP98, we discussed a method to obtain informations on the cosmological 
parameters £l m and £l\ while reconstructing the lens gravitational potential of clusters with 
multiple image systems at different redshifts. 
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Figure 3: Left. Generation of images by a zi = 0.3 cluster with ao = 1400 km/s, # a = 13.54" and 9 S — 145.8". 
Close to their respective critic lines, we see 3 families of images at z s i = 0.6, z S 2 = 1 and z S 3 = 2. Right. Solid lines 
: X 2 (^rn, confidence levels obtained for this configuration. Generating arbitrary values: (fi^n, ) = (0.3, 0.7). 
Dashed lines: x 2 (^m, fix) confidence levels obtained considering 10 clusters in this same configuration. 




Figure 4: HST image of MS2137-23 (zi = 0.313). We see 3 families of images, a tangential arc A1-A2 with its 
associated arcs A3 and A4, and 2 radial arcs Bl and CI with their associated images B2 and C2. 



This technique gives degenerate constraints, fl m and £1\ being negatively correlated, with a 
better constraint of the matter density. With a single cluster in a typical lensing configuration 
we can expect the following error bars : Vt m = 0.3±0.24, £l\ = 0.7±0.5. To perform that, several 
general conditions must be fulfilled: 

* a cluster with a rather regular morphology, so that a few parameters are needed to describe 
the gravitational potential ; this is not so restrictive because we saw that a bimodal cluster can 
also provide a constraint, 

* "numerous" systems of multiple images, probing each part of the cluster, 

* a good spatial resolution image (HST observations) of the cluster and arcs to compute 
relatively precise - 0.1" - (flux weighted) images positions, 

* spectroscopic precision for the different redshifts that should be also regularly distributed 
from %\ to high values - this requires deep spectroscopy on 8-10m class telescopes due to the 
faintness of the multiple images . 

It is important to notice that one cluster could provide one constraint on the geometry of 
the whole universe. And it is possible to combine data from different clusters to tighten the 
error bars. Combining the study of about 10 clusters would lead to meaningful constraints. The 
dashed lines confidence levels in the Fig. || are the result of a numerical simulation made with 
10 identical clusters. 

Actually the degeneracy depends only on the different redshifts involved that we will have 
various sets of when applying the method to real configurations. This should lead to a more 
reduced area of allowed cosmological parameters. We are encouraged by more and more known 
observations including systems with multiple sources and we plan to apply this technique to 
clusters like MS2137-23, MS0440+02, A370, AC114 and A1689. 
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